
IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

(Docket No. 15-XZ-4971) 




PATENT APPLICATION OF: 



KENNETH KUMP 



EXPRESS MAIL NO. FJ.282688023US 



SERIAL NO.: 



Date: June 24. 1999 



FILED: 

FOR: METHOD AND APPARATUS FOR 

DETERMINING A DYNAMIC RANGE 
OF A DIGITAL MEDICAL IMAGE 



Assistant Commissioner for Patents 
BOX PATENT APPLICATION 
Washington, D.C. 20231 

Dear Sirs: 

Transmitted herewith for filing is the patent application of Inventor (s): Kenneth Kump. 



[X] A Combined Declaration and Power of Attorney signed by each inventor. 
[X] 5 sheet(s) of drawing(s) (informal). 

[X] 23 pages of specification including a 1 page abstract, and 5 page(s) of claims 
(Claims 1-24). 

[ ] A verified statement to establish small entity status under 37 C.F.R. 1 .9 and 37 
C.F.R. 1.27. 

[ ] An associate power of attorney. 

[ ] A Transmittal of Information Disclosure Statement and its enclosures. 
[] Priority Document 



Enclosed are: 



2. Benefit of the filing date of each priority document or U.S. provisional patent application 
hsted below (if any) is claimed under 35 U.S.C. Section 1 19 and Section 120: 

Appl'n. 

Docket No. Country Serial No. Filing Date 



3, The filing fee has been calculated as shown below, based on the assignee's status as a 
large entity: 

Claims Extra Rate per 

Presented Claims Extra claim 



Total Claims 24 -20 4 X $ 1 8 $72.00 

Indep Claims 2-3 0 X$72 $.00 

Basic fee $760.00 

Multiple dependent claim presented ($260 if any present) £00 

Total Filing fee $^32.00 



4. The following arrangements have been made to pay the filing fee: 

[X] A check in the amount of $832.00 to cover the filing fee is enclosed. 

[ ] The filing fee is not enclosed. This form does not authorize a charge to a deposit 
account. 

[ ] Please charge our Deposit Account No. 13-0017 in the amount of $--. 

[X] The Commissioner is hereby authorized to charge any additional filing fees which 
may be required under 3 7 C.F.R. 1 . 1 6 for this application and for any amendment 
thereto, or credit any overpayment to Account No. 13-0017. 

5. The Commissioner is hereby authorized to charge payment of the following fees during 
the pendency of this appUcation or credit any overpayment to Deposit Account No. 13-0017. 

[ ] No charges are authorized. 

[X] Any patent appUcation processing fees under 37 C.F.R. 1.17. 

6. Please address all telephone calls to: 

at telephone No. (312) 707-8889. 



7. 



Please address all correspondence to: 



McAndrews, Held & Malloy, Ltd. 
34th Floor 
500 West Madison Street 
Chicago, Illinois 60661 



X3 J «^50W~-- 

nate: .Tune 24. 1999 



Dean D. Small 
Reg. No. 34,730 



C3 
IJ 

□ 

i-i 

a 
m 
ru 
+ 



-3- 



UNITED STATES PATENT APPLICATION FOR 

METHOD .AND APPARATUS FOR DETERMINING A 
DYNAMIC RANGE OF A DIGITAL MEDICAL IMAGE 

Inventor (s) : 

Kenneth Kump 



Prepared by: 
McAndrews, Held & Malloy 
500 West Madison Street 

Suite 3400 
Chicago, Illinois 60661 

(312) 707-8889 



Docket No. 15-XZ-4971 



EXPRESS MAIL 

^^Express Mail" mailing label number EL282688Q23US 
Date of Deposit: June 24, 1999 



TITLE OF THE INVENTION 

METHOD AND APPARATUS FOR DETERMINING A 
DYNAMIC RANGE OF A DIGITAL MEDICAL IMAGE 

CROSS-REFERENCE TO RELATED APPLICATIONS 
Not applicable. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 

DEVELOPMENT 

Not applicable. 

BACKGROUND OF THE INVENTION 
The present invention generally relates to medical 
diagnostic imaging and in particular to a method and 
apparatus for determining dynamic range of a digital 
medical image to be displayed. 
5 X-ray imaging had long been an accepted medical 

diagnostic tool. X-ray imaging systems are commonly used 
to capture, as examples, thoracic, cervical, spinal, 
cranial and abdominal images that often include the 
information necessary for a doctor to make an accurate 

10 diagnosis. When having a thoracic x-ray image taken, for 
example, a patient stands with his or her chest against 
an x-ray sensor as an x-ray technologist positions the x- 
ray sensor and an x-ray source at an appropriate height. 
The x-ray sensor then detects the x-ray energy generated 

15 by the source and attenuated to various degrees by 
different parts of the body. An associated control 
system scans out the detected x-ray energy and prepares a 
corresponding diagnostic image on a display. Optionally, 
the x-ray sensor may be a solid state digital image 
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detector. If the x-ray sensor is a conventional 

screen/film configuration, the screen converts the x-rays 
to light, to which the film is exposed. 

In conventional radiographic imaging systems, the x- 
5 ray technique is chosen by the operator. The operator or 
the automatic exposure control system selects or 
determines a desired exposure for the selected 
screen/film configuration in order to obtain a desired 
optical density of the exposed film. The optical density 

10 represents the "lightness" or "darkness'^ of the screen, 
detector or film once exposed to x-rays . By controlling 
the manner (e.g., time, orientation, etc.) of exposure by 
the detector, screen or film to x-rays, the film 
lightness or darkness may be varied. It is preferable to 

15 achieve a consistent optical density from one exposure to 
the next in order to facilitate diagnosis and examination 
by physicians when analyzing radiographic images. 
Different exposures arise from one patient to the next, 
from one film type to the next, from one medical imaging 

2 0 system to the next, from one orientation to the next and 

the like. 

In the past, it has been quite difficult to maintain 
a uniform optical density from one exposure to the next 
(e.g., patient to patient, film to film, system to 
25 system, patient angle to patient angle) due to inherent 
differences. For instance, each patient has a slightly 
different size and anatomy which causes the internal 
organs of the patient to be located at different 
positions relative to the detector or screen/film. For 

3 0 example, when attempting to obtain an x-ray of a chest 

image, every patient's lungs and rib cage are of a 
different size. The position of the lungs is also 
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somewhat unknown which creates a large variance in the 
resulting exposure. Further, patient position is not 
precisely controlled and hence each patient is located in 
a slightly different position or orientation with respect 
5 to the detector or screen/film configuration. Variation 
in patient position and orientation further create 
variance in the resulting exposure. Optical density may 
further be varied due to the particular pathology 
followed by the x-rays through the patient, due to 

10 foreign objects within a patient (e.g., pacemakers and 

the like) as well as due to differences in patient 
thickness and resulting scatter pattern properties. 

An automatic exposure control has been proposed for 
use with radiographic systems in an attempt to control 

15 the optical density of the exposed film. Automatic 
exposure control systems typically use an x-ray sensitive 
ion chamber located proximate the detector, screen/film 
configuration and arranged to be proximate a particular 
anatomy of a patients during examination . For instance , 

2 0 an ion chamber may be located within a region of the 
detector or screen/film configuration calculated to be 
proximate the patient's lung during a particular form of 
examination. Alternatively, or in addition, an ion 
chamber may be located proximate the patient's 

2 5 mediastinum. The automatic exposure control measures the 

x-rays detected by the ion chamber and terminates the 
exposure when a preset dose is measured. 

However, automatic exposure control systems have 
experienced difficulties. In particular, the position of 

3 0 an individual patient's lung is unknown at the time that 

the ion chamber is placed proximate the detector, screen 
or film. Hence, different patients continue to create a 
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large variance in the resulting exposure to the ion 
chamber. For instance, the ion chambers may not actually 
be located proximate certain patient's lungs or 
mediastinum. When an ion chamber is located proximate an 
5 anatomy other than the lung or mediastinum, the automatic 
exposure control terminates exposure based on inaccurate 
measurements. A certain percentage of chest films result 
in creation of either too dark or too light of an image. 
When the image is too dark or too light, it may be 

10 necessary to repeat the x-ray examination to retake the 
medical image. It is quite time consuming to retake 
medical images. Film development may require a 

relatively long period of time, such as five to fifteen 
minutes, during which the patient may leave the image 

15 acquisition area. 

Further, a resulting presentation of a medical image 
is determined by the selection of the type of detector, 
film/screen configuration in combination with the desired 
x-ray technique. Different types of detectors and 

20 screens/film configurations experience different amounts 
of image noise. In the past, noise has been partially 
corrected by varying the input exposure time. However, 
to maintain a constant optical density from one exposure 
to the next, when detector, film or screen types are 

2 5 changed, the exposure time must be changed in order to 

account for the fixed dynamic range of the new detector, 
screen/film configuration. It is quite cumbersome to 
change detectors, screens or films, and thus rarely done. 

More recently, digital detectors have been proposed 

3 0 for use with radiographic imaging. Digital detectors 

afford a significantly greater dynamic range than 
convention screen/film configurations, typically as much 
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as two to three times greater. Heretofore, the automatic 
exposure control and operator must still be relied upon 
to limit the exposure of the digital detector to account 
for the detectors greater dynamic range. 

A need remains for an improved dynamic range 
detection and control method and apparatus for use with 
digital medical imaging, such as in radiographic imaging. 

SUMMARY OF THE INVENTION 

A preferred embodiment of the present invention 
provides a method and apparatus for determining a dynamic 
range of a digital medical image for a medical imaging 
system. The medical image contains clinical regions and 
non-clinical regions. The method and apparatus identify 
the non-clinical regions of the digital medical image and 
mask the non-clinical regions from the digital medical 
image to provide an automatic contrast control for 
displaying a clinical image. Once the non-clinical 
region has been masked, the method and apparatus 
determine the dynamic range for the clinical image. 

According to one alternative embodiment, the non- 
clinical regions are identified by dividing the digital 
medical image into bands of predetermined width and 
locating the non-clinical region in each band. The bands 
may be aligned horizontally and/or vertically within the 
digital medical image. During identification, the bands 
of the digital medical image are differentiated and 
discriminated with thresholds to identify changes in gray 
scale along a profile corresponding to a horizontal or 
vertical band. A mask for masking non-clinical regions 
is turned on and off based on maximum and minimum peaks 
in the differentiated profile for the digital medical 
image . 
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In an alternative eTnbodiment , the non- clinical 
regions are identified after a histogram is generated 
from the digital medical image. Thresholds may be 
obtained based on a dynamic range of the digital medical 
5 image and used to identify maximum and minimum gray scale 
values in the histogram associated with non-clinical 
regions. The non-clinical regions are masked based on 
the maximum and minimum values identified in the 
histogram as being associated with non-clinical regions. 

10 In yet a further alternative embodiment of the 

present invention, a medical diagnostic imaging system is 
provided for determining a dynamic range of a digital 
medical image to be displayed. A digital detector may be 
provided in order to obtain the digital medical image. 

15 The digital medical image includes a clinical region and 
a non-clinical region. The system includes a 

segmentation module and a dynamic range module. The 
segmentation module identifies the non-clinical region in 
the digital medical image . The dynamic range module 

2 0 determines the dynamic range of the clinical region of 
the digital medical image once the non-clinical region 
has been segmented. The segmentation module may identify 
raw radiation regions and/or collimated regions as non- 
clinical regions. 

25 In one alternative embodiment, the segmentation 

module identifies non-clinical regions based on 
variations in gray scale levels of the digital medical 
image. The segmentation module differentiates and 

threshold detects at least a portion of the digital 

30 medical image to identify the non-clinical region. In an 
alternative embodiment, the segmentation module generates 
a histogram corresponding to the digital medical image 
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and utilizes at least one gray scale threshold value to 
discriminate the non-clinical region. 

A processor may be included to calculate at least 
one threshold based on a dynamic range of the digital 
5 medical image. The dynamic range module may include a 
processor masking non-clinical regions when determining 
the dynamic range of the clinical region. A processor 
may be further provided to calculate at least one of a 
maximum and minimum gray scale level for the digital 

10 medical image in order to identify the non-clinical 
region. A processor may also calculate at least one of a 
maximum and minimum gray scale level for the clinical 
region in order to determine its dynamic range. 

At least one preferred embodiment of the present 

15 invention frees the reliance of the system upon patient 
positioning and x-ray acquisition techniques in order to 
obtain a uniform output density for display from patient 
to patient, film to film and system to system. The 
dynamic range management methods and systems of the 

2 0 preferred embodiments of the present invention are 
equally useful in general radiographic applications and 
may be scalable with respect to various exposures and 
mean gray scale levels. A resulting image presentation 
may improve signal-to-noise ratio characteristics without 

2 5 changing a global contrast of an image, thereby taking 

advantage of the wide dynamic range of digital detectors. 

By affording the ability to detect the 
diagnostic/clinical dynamic range of a digital medical 
image, specific knowledge is no longer required of the 

3 0 dose level, nor patient type. Optionally, dose level 

and/or patient type may be used to improve robustness. 
The dynamic range management methods and systems of the 
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preferred embodiments of the present invention allow 
independent operation of image presentation and image 
acquisition, while assuring image quality and consistent 
presentation from image to image, patient to patient, 
dose to dose and system to system. The preferred 
embodiments of the present invention also afford the 
ability to simulate global contrast curves of screen/ 
films . 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates a preferred embodiment of the 
present invention . 

Fig. 2 illustrates a processing sequence carried out 
according to a preferred embodiment of the present 
invention. 

Fig. 3 illustrates a digital medical image divided 
in vertical bands generated in accordance with a 
preferred embodiment of the present invention. 

Fig. 4 illustrates a digital medical image divided 
into horizontal bands in accordance with a preferred 
embodiment of the present invention. 

Fig. 5 illustrates a gray scale profile of a 
horizontal or vertical band of a digital medical image 
generated in accordance with a preferred embodiment of 
the present invention. 

Fig. 6 illustrates a differentiated gray scale 
profile associated with the profile of Fig. 5 generated 
in accordance with a preferred embodiment of the present 
invention. 

Fig. 7 illustrates a processing sequence carried out 
in accordance with an alternative embodiment of the 
present invention. 
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Fig. 8 illustrates a gray scale histogram generated 
in accordance with an alternative embodiment of the 
present invention . 

DETAILED DESCRIPTION OF THE INVENTION 

Fig. 1 illustrates a medical imaging system 
5 configured in accordance with the preferred embodiment of 
the present invention. The medical system 8 includes a 
dynamic range determining system 10 and a patient 
examining subsystem 12 . The patient examining subsystem 
12 includes a source 16 controlled by a control circuit 

10 14. The source 16 emits a medium, such as radiation, 
through a patient 20. A digital detector 22 detects the 
medium 18 that passes through the patient 20. A detector 
output circuit 24 converts the information stored by the 
digital detector 22 to an electronic format processable 

15 by the processor 2 6 in the dynamic range control 
determining 10. The processor 26 accepts an input from 
the detector output circuit 24 and based thereon stores a 
raw digital image 28 in memory 27. The raw digital image 
28 includes a clinical region and a non-clinical region. 

20 The processor 26 distinguishes the non-clinical 

region from the clinical region in the raw image 2 8 in 
accordance with alternative techniques described below to 
form a segmented image 30. The segmented image 30 is 
stored in memory 27 and contains only the clinical 

25 region, with the non- clinical region being masked. Next, 
the processor 2 6 calculates the dynamic range of the 
clinical region in the segmented image 3 0 and based 
thereon generates dynamic range image characteristics 32. 
The dynamic range image characteristics 32 are stored in 

3 0 memory 2 6 . 
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Turning to Fig. 2, a processing sequence is 
illustrated that may be carried out by the processor 2 6 
in order to determine the dynamic range of the digital 
medical image. Beginning at step 50, the processor 2 6 
5 obtains the raw digital image 2 8 either from detector 24 
or memory 27, At step 52, the processor 26 divides the 
raw digital image 28 into horizontal and/or vertical 
bands of pixels. Exemplary bands are illustrated in 
Figs. 3 and 4. Fig. 3 illustrates a raw digital image 7 0 

10 which may be comprised of an array of pixel values, such 
as a 2K by 2K pixel array. The digital image 70 is 
divided at step 52 into a plurality of vertical bands 72 
extending from the top 73 to the bottom 75 of the image 
and having a predetermined vertical bandwidth 74. By way 

15 of example only, each vertical bandwidth 74 may equal an 
even number of pixels, such as 100 in vertical band rows 
83 -85 . Thus , in a 2000 by 2 000 digital image 70 , twenty 
vertical bands 72 would be utilized, each of which is 100 
pixels in width. 

20 Alternatively or in addition, at step 52, the 

digital image 70 may be divided into horizontal bands 76 
as illustrated in Fig. 4. The horizontal bands 76 extend 
from the left side 77 to the right side 79 of the digital 
image 76 and have a predetermined horizontal bandwidth 

25 78. By way of example, the horizontal bandwidth 78 may 
also be 100 pixels in width horizontal band columns 80- 
82 . The horizontal bands 76 extend from one side of the 
digital image 70 to the other. 

Once the horizontal and/or vertical bands 72, 76 are 

3 0 created in step 52, the processor 26 calculates an 
average gray scale pixel value across the width of each 
horizontal and vertical band 72 and 76 (step 54) . By way 
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of example, a horizontal band 78 may be 100 pixels in 
width and 2, 000 pixels in length. At step 54, the 
processor sequentially steps through each set of 100 
pixels in horizontal band columns 80-82. An average 
5 pixel value is calculated for the 100 pixels in 
horizontal band column 80, while separate average pixel 
values are calculated based upon the 100 pixel sets in 
horizontal band columns 81 and 82. This process is 
repeated until a subsampled vector is formed for each 

10 horizontal band 78. In the example of Fig. 4, each sub- 
sampled vector includes 2 , 000 average pixel values 
extending across the profile or length of the horizontal 
band 78. A separate sub-sampled vector is calculated for 
each horizontal band 78 and/or each vertical band 74. 

15 Once the sub- sampled vectors are calculated, a horizontal 
and/or vertical profile is generated for each horizontal 
and vertical band 74 and 78. 

Fig. 5 illustrates graphically a profile 90 
associated with an exemplary vertical or horizontal band 

20 74 or 78. The horizontal axis 92 of the profile 90 in 

Fig. 5 corresponds to the length of an associated 
vertical or horizontal band 74 or 78 . The vertical axis 
94 corresponds to the gray scale levels of the average 
pixel values for a subsample vector associated with a 

2 5 grouping of pixels within one of rows 83-85 or columns 
80-82 of an associated vertical or horizontal band 74 or 
78, respectively. In the example discussed above, where 
a horizontal band 78 includes 2,000 pixels in length and 
100 pixels in w^idth, the associated profile 90 would 

30 include 2,000 average pixel values spanning the length of 
horizontal axis 92 of the profile 90. 
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The exemplary profile 90 of average pixel values 
illustrated in Fig. 5, may be from an x-ray scan of the 
chest of a patient . The profile 90 includes leading and 
trailing portions 95 and 110 having gray scale levels 
5 approaching zero. Leading and trailing portions 95 and 
110 may correspond to regions of the detector covered by 
a collimator which blocks a majority of the radiation. 
Raw radiation peaks 96 and 108 correspond to areas having 
high exposure to radiation and may simply represent raw 

10 radiation exposed to the digital detector. Tissue 
valleys 98 and 106 correspond to tissue portions of the 
patient such as along the left and right sides of the 
chest. Intermediate lung peaks 100 and 104 correspond to 
the lungs of the patient, while mediastinum valley 102 

15 corresponds to the center or mediastinum of the patient. 

Once the profile 90 has been generated, the profile 
90 of the sub- sample vector for each band is 
differentiated at step 56 to obtain a differentiated gray 
scale profile identifying the change in gray scale at 

20 each location along profile 90. Fig. 6 illustrates a 
differentiated gray scale profile 120 with the horizontal 
axis 122 corresponding to a length of an associated 
vertical or horizontal band 74 or 78 . The vertical axis 
124 of the differentiated gray scale profile 120 

25 identifies positive and negative changes in the gray 
scale level per unit distance along an associated band 
72, 74. For example, peak 12 6 is associated with the 
point 97 of profile 90 having a maximum slope. The peak 
99 in the raw radiation portion 96 corresponds to point 

30 128 in the differentiated gray scale profile 120, where 

the differentiated gray scale equals zero. The points 
130-136 at which the differentiated gray scale profile 
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120 cross zero correspond to points 140-146 in the gray 
scale profile 90, respectively. 

Next, at step 58 in Fig. 2, the processor 26 
calculates thresholds that may be utilized to identify 
5 peaks in the differentiated gray scale profile 120. At 
step 60 , the processor 2 6 discriminates each 
differentiated band to identify starting and stopping 
points of non-clinical regions, such as associated with 
raw radiation and collimators. By way of example only, a 

10 raw radiation search algorithm may be carried out to 
analyze the differentiated gray scale profile 120 to 
identify the starting and stopping points of the non- 
clinical regions . 

As an example, the processor 2 6 may begin by 

15 scanning the differentiated gray scale profile 12 0 from 
left to right and identify the first positive peak 12 6 of 
the profile 120. Once peak 126 is identified, the 
processor 2 6 reverses scanning direction and searches 
toward the left edge of the profile 12 0 until the 

20 processor 26 identifies a point where differentiated 
vector equals zero. The processor 26 then turns on a 
mask in the digital image 2 8 at the point corresponding 
to the point where the differentiated vector equals zero 
(namely point 12 5) . Next, the mask remains on as the 

25 differentiated profile is scanned for the first negative 
peak (point 12 9) . Peaks 12 6 and 12 9 may be identified by 
comparing the peak thresholds 121 and 123 to the 
differentiated profile 120. Once negative peak 129 is 
identified, the processor 2 6 continues scanning the 

3 0 differentiated vector until it equals zero (namely point 
13 0) - The mask is turned off at point 13 0. 
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The processor 12 6 continues to scan the 
differentiated profile 120 until the profile 120 again 
exceeds the positive peak threshold 121 (such as at point 
150) . When peak 150 is identified, the processor 26 
5 reverses scanning direction to search until it identifies 
the point where the differentiated vector equals zero 
(point 134) and again turns on the mask. Scanning again 
continues to the right while the processor 26 searches 
for a negative peak which exceeds the threshold 123 (such 

10 as at point 151) . Once negative peak 151 is identified, 
the differentiated vector 12 0 is then scanned for the 
next point at which it equals zero (point 136) . At point 
13 6, the mask is again turned off. 

The foregoing discrimination technique effectively 

15 removes the portion of the digital medical image 
preceding point 130 and following point 134. 

Turning to Fig. 5, the above discussed 
discrimination technique has the effect of removing or 
ignoring the portion of the digital medical image 

20 preceding point 140 and following point 144. The 
processor effectively identifies everything preceding 
point 14 0 and following point 144 to correspond to non- 
clinical regions, namely raw radiation and collimator 
regions. The remaining section of the profile 90 between 

25 points 140 and 144 include lung peaks 100 and 104 which 
correspond to the lungs and a mediastinum valley 102 
which corresponds to the mediastinum. The region of the 
profile 90 between points 140 and 144 corresponds to the 
segmented image 3 0 stored in memory 2 7 in the dynamic 

30 range determining system 10. 

Next, at step 64, the processor 2 6 computes desired 
image characteristics, such as mean, median, average. 
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Standard deviation, maximum and minimum gray scale values 
of the segmented image 30 and the like. As shown in Fig. 
5, the maximum and minimum gray scale values correspond 
to points 141-143 and define the dynamic range of the 
clinical region. After determining the desired image 
characteristics, such as maximum and minimum gray scale 
values of the segmented image 30, in step 66, the 
processor 26 adjusts the dynamic range of the medical 
diagnostic imaging system to form a dynamic range 
adjusted image based on the image characteristics. The 
dynamic range adjusted image and/or dynamic range image 
characteristics 32 are stored in memory 26. 

Optionally, the dynamic range adjusted image may be 
obtained by mapping the clinical region of the original 
raw digital image 28 into an output image segment by 
passing the clinical region of the raw digital image 2 8 
through a look-up table that correlates input pixel 
values to desired output pixel values. For instance, the 
input and output pixel values stored in the look-up table 
may have a linear relation to one another, a sigmoidal 
relation to one another and the like. The image 

characteristics, such as maximum and minimum gray scale 
values, calculated at step 64 may be utilized to adjust 
the slope or shape of the mapping function in the look-up 
table, in order to maintain a constant dynamic range for 
all digital images processed by the system regardless of 
variations in film, patients, exposure and systems. 

It should be noted that not all digital medical 
images include non-clinical regions. Hence, at step 60, 
the discrimination may have a null or empty solution 
(e.g., no starting and stopping points). Digital medical 
images that do not include non-clinical regions will be 
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treated in step 64 as the segmented image. Step 62 will 
effectively be skipped. 

Next , the discussion turns to an alternative 
embodiment for segmenting non-clinical regions, as 
5 illustrated in Fig. 7. Beginning at step 200, the 
processor 2 6 again obtains a raw digital image 28. A 
histogram is generated at step 2 02 from the raw digital 
image 28. At step 204, gray scale thresholds are 
calculated for gray scale levels within the histogram 

10 associated with non-clinical regions. At step 206, the 
gray scale thresholds are used to identify peaks in the 
histogram which correspond to the non-clinical regions. 
Once the non-clinical regions are identified in step 206, 
the gray scale levels of the histogram that are 

15 associated with the non-clinical are masked and the 
remaining gray scale levels from the digital image form a 
segmented image 3 0 corresponding to the clinical region. 
At step 210, the processor determines dynamic range 
characteristics for the segmented image 30. At step 212, 

2 0 the dynamic range of the medical diagnostic imaging 

system is adjusted to form a diagnostic range adjusted 
image based on the dynamic range characteristics 
calculated in step 210. 

The processing sequence of Fig. 7 is better 
25 illustrated in connection with Fig. 8. Fig. 8 

illustrates a histogram 300 generated at step 202, in 
which the gray scale levels are denoted along the 
horizontal axis 3 02 and the count of pixels at each gray 
scale level is denoted along the vertical axis 3 04. The 

3 0 histogram 3 00 includes a leading region 3 06 having gray 

scale levels that are very low. Portion region 3 06 may 
correspond to a collimator region. The histogram 3 00 
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further includes a trailing region 308 having a large 
number of pixels at very high gray scale levels. Region 
3 08 may correspond to raw radiation regions of the 
digital image 28. An intermediate region 310 corresponds 
5 to the clinical region containing the gray scale levels 
associated with the lungs and mediastinum. 

During operation, gray scale thresholds 312 and 314 
are calculated at step 204 and used to identify peaks 
corresponding to non-clinical regions (i.e., regions 306 

10 and 308) . At step 208, the gray scale levels in regions 
306 and 308 are masked to hide associated non-clinical 
regions, and to form a segmented image 3 0 containing only 
clinical regions have gray scale values from intermediate 
region 310. Next, at step 210, the maximum and minimum 

15 gray scale values within intermediate region 310 are 
calculated and used to adjust the dynamic range of the 
medical system at step 212. 

The dynamic range management system of the preferred 
embodiments enables desired output optical densities to 

2 0 be obtained for all patients, regardless of the type of 

film, screen or detector used, regardless of the 
exposure, regardless of the dose and regardless of the 
individual patient being scanned. 

While particular elements, embodiments and 
25 applications of the present invention have been shown and 
described, it will be understood, of course, that the 
invention is not limited thereto since modifications may 
be made by those skilled in the art, particularly in light 
of the foregoing teachings. It is therefore contemplated 

3 0 by the appended claims to cover such modifications as 

incorporate those features which come within the spirit 
and scope of the invention. 
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CLAIMS : 

1 . A method for determining a dynamic range of a 
digital medical image for a medical imaging system, the 
digital medical image containing a clinical region, 
comprising : 

determining whether the digital medical image 
includes a non-clinical region; and 

calculating a dynamic range for a clinical region. 

2. The method of claim 1, further comprising: 
dividing the digital medical image into bands of 

predetermined width, said determining step identifying 
the non-clinical region in each band. 

3. The method of claim 1, wherein the digital 
medical image is divided into horizontal and vertical 
bands during the determining step. 

4. The method of claim 1, further comprising: 
differentiating said digital medical image, said 

determining step calculating a position of the non- 
clinical region based on a result of said 
differentiation. 

5. The method of claim 1, further comprising: 
calculating at least one threshold based on a 

dynamic range of the digital medical image, said at least 
one threshold being used to identify at least one of 
maximum and minimum values for the non-clinical region. 

6. The method of claim 1, wherein said identifying 
step discriminates at least one of histogram maximum and 
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minimum values for a non-clinical region based on at 
least one predetermined threshold. 

7. The method of claim 1, further comprising: 
when a non-clinical region is determined to exist, 

5 masking the non-clinical region from the digital medical 
image before calculating said dynamic range. 

8. The method of claim 1, further comprising 
masking non-clinical regions based on at least one of 
gray scale maximum and minimum values for the non- 
10 clinical region. 

9. The method of claim 7, further comprising: 
generating a histogram of the digital medical image, 

said masking step masking gray scale levels from the 
histogram that exceed predetermined upper and lower 
15 thresholds. 

10. The method of claim 1, wherein said determining 
step determines that the digital medical image does not 
include a non-clinical region and said calculating step 
calculates a dynamic range of the entire digital medical 

2 0 image as the clinical region . 

11. A medical diagnostic imaging system for 
controlling a dynamic range of a digital medical image to 
be displayed, the digital medical image including a 
clinical region and a non-clinical region, comprising: 

25 a segmentation module identifying a non-clinical 

region of a digital medical image; and 
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a dynamic range module determining a dynamic range 
of a clinical region of the digital medical image once 
the non-clinical region is segmented. 

12. The system of claim 11, further comprising: 

5 a digital detector obtaining a digital medical image 

having clinical and non-clinical regions. 

13. The system of claim 11, wherein the 
segmentation module identifies a raw radiation region in 
the non-clinical region. 

10 14. The system of claim 11, wherein the 

segmentation module identifies a collimated region in the 
non-clinical region. 

15. The system of claim 11, wherein the 
segmentation module identifies non-clinical regions based 

15 on variations in gray scale levels of the digital medical 
image . 

16. The system of claim 11, wherein the 
segmentation module differentiates at least a portion of 
the digital medical image to identify the non-clinical 

2 0 region. 

17. The system of claim 11, wherein the 
segmentation module discriminates the non-clinical region 
based on at least one gray scale threshold value. 

18. The system of claim 11, further comprising a 
25 processor calculating at least one threshold based on a 
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dynamic range of the digital medical image, said 

segmentation module discriminating the non-clinical 
region based on said threshold. 

19. The system of claim 11, wherein said dynamic 
5 range module including a processor masking over a non- 
clinical region when determining the dynamic range of the 
clinical region. 

20. The system of claim 11, further comprising: 

a processor calculating at least one of a maximum 
10 and minimum gray scale level for the digital medical 
image in order to identify the non- clinical region. 

21. The system of claim 11, further comprising: 

a processor calculating at lest one of maximum and 
minimum gray scale levels for the clinical region in 
15 order to determine the dynamic range of the clinical 
region. 

22. The system of claim 11, further comprising: 

a processor generating a histogram of at least a 
portion of the digital medical image to identify gray 
2 0 scale levels associated with non-clinical regions. 

23. The system of claim 11, wherein the 
segmentation module masks a non-clinical region 
identified in the digital medical image. 

24. The system of claim 11, wherein the 
25 segmentation module determines that the digital medical 

image does not include a non-clinical region, said 
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dynamic range module determining a dynamic range of the 
digital medical image as the clinical region. 
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METHOD AND APPARATUS FOR DETERMINING A 
DYNAMIC RANGE OF A DIGITAL MEDICAL IMAGE 

ABSTRACT OF THE DISCLOSURE 

A method and apparatus are provided for adjusting a 
dynamic range of a digital medical image for a medical 
imaging system. The digital medical image contains a 
clinical region and a non- clinical region . The method 
5 and apparatus identify the non-clinical region of the 
digital medical image and mask the non-clinical region 
therefrom to form a clinical image. The clinical image 
is then used to calculate a desired dynamic range for the 
medical imaging system. The dynamic range of the medical 

10 imaging system is adjusted accordingly. According to one 
embodiment, the non-clinical region is identified by 
dividing the digital medical image into bands of a 
predetermined width, generating profiles for each band 
and differentiating the profiles to obtain a 

15 differentiated profile of each band of a digital medical 
image. The differentiated profile is than analyzed to 
identify peaks that exceed predetermined thresholds, 
wherein the regions of the differentiated profile 
proximate the peaks exceeding the threshold correspond to 

20 non-clinical regions. Once the non-clinical regions are 
identified, they are masked or removed . Next , a desired 
image characteristic, such as maximum and minimum gray 
scale values, are determined for the clinical region and 
a desired dynamic range for the image is obtained based 

25 on the image characteristics of the clinical region. In 
an alternative embodiment, a histogram is used to 
identify the non-clinical regions which are subsequently 
masked from the digital medical image. 
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I believe I am the origiml, first and sole inventor (if only one name is listed below) or an 
original, first and joint inventor (if plural names are listed below) of the subject matter which is 
claimed and for which a patent is sought on the invention entitled 

METHOD AND APPARATUS FOR DETERMINING A DYNAMIC RANGE OF A 

DIGITAL MEDICAL IMAGE 

the specification of which: 

X is attached hereto. 

was filed on - as United States Application Number or PCX Intemational 
AppUcation Number ^ and was amended on 

I hereby state that I have reviewed and understand the contents of the above identified 
specification, including the claims, as amended by any amendment referred to above. 

I acknowledge the duty to disclose to the United States Patent and Trademark Office all 
infoimation which I know to be material to patentability as defined in Title 37, Code of Federal 
Regulations, §1.56. 

I hereby claim foreign priority benefits under Title 35, United States Code, §119(a)-(d) of any 
foreign application(s) for patent or inventor^s certificate listed below, I have also identified 
below any foreign application for patent or inventor's certificate having a filing date before that 
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